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Ashok Das 

Department of Physics and Astronomy, University of Rochester, Rochester, NY, USA 

J. Frcnkel and S. H. Pereira 
Instituto de Fisica, Universidade de Sao Paulo, Sao Paulo, Brazil 

J. C. Taylor 

DAMTP, Centre for Mathematical Sciences, Cambridge University, Cambridge, UK 

We study some aspects of the quantum theory of a charged particle moving in a time-independent, 
uni-directional magnetic field. When the field is uniform, we make a few clarifying remarks on the 
use of angular momentum eigenstates and momentum eigenstates with the diamagnetism of a free 
electron gas as an example. When the field is non-uniform but weakly varying, we discuss both 
perturbative and non-perturbative methods for studying a quantum mechanical system. As an 
application, we derive the quantized energy levels of a charged particle in a Helmholtz coil, which 
, go over to the usual Landau levels in the limit of a uniform field. 
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^ 1 I. INTRODUCTION 

> \ 

^>0 \ Although the quantum theory of a charged particle moving in a uni-directional constant magnetic field has 
been studied over many decades (for some text-book treatments, see for example pj), we would like to make 
some clarifying remarks which we believe are useful. In studying this problem, one has to make a choice 
of the gauge. Different gauge choices, of course, lead to a trivial phase redefinition of the wave function. 
Moreover, corresponding to specific gauges, one chooses either a wave function which is an eigenstate of the 
momentum operator (which has continuous eigenvalues) or that of the angular momentum operator (which has 
discrete eigenvalues) . In section II we point out that such a choice is unessential and derive an explicit relation 
between the two bases. This is then used in the treatement of diamagnetism of a free electron gas with angular 
momentum eigenstates. In section III, we study the case when the magnetic field has a weak z-dependence 
(where the z-axis is in the direction of the field). Then, there must be also a weak radial field and we work 
under conditions where it is consistent to treat this component of the field perturbatively. We find the energy 
eigenvalues for certain models of the z-dependence. These models could approximate the field between two 
parallel circular currents. In section IV, we discuss the form of the wave-function for large quantum numbers 
associated with the z-motion, using a non-perturbative approach in a slowly varying field. As an application, 
we study in section V the structure of the quantized energy levels of a charged particle in a Helmholtz coil, 
which are specified by two quantum numbers. The first is associated with the motion in the (x,y) plane, 
while the second characterizes the motion in the z-direction. This energy spectrum may be regarded as a 
generalization of the Landau levels [l| to the case of a slowly varying magnetic field. A brief conclusion is 
presented in section VI. 

II. UNIFORM MAGNETIC FIELD 

We take the field to be in the z-direction, and to have z-component B. We take the particle's charge and 
mass to be e and to. Two gauges for the vector potential which are commonly used are 

A = {—By/2, Bx/2, 0) (1) 

and 

A = (—By, 0, 0). (2) 

But the difference between these two gauges can only be trivial, since the wave- functions just differ by the 
phase factor 

u — exp[ieBxy / (2ch)]. (3) 
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Therefore we might just as well choose {Q. 
For shortness, we define 



B = eB/c 



The Hamiltonian in gauge (JIJ is 



H= 2m~ [{Px ~ BV,2) + {Py + Bx/2) + P * ] 



(4) 



where 



1 

2m 



£- [Pl + vl +Pl+ B 2 (x 2 + y 2 )/4 + BL] = H 1 + p 2 z /2m, 



k d k d -t, d 

Px = -iri—, Py = -in—, p z = -in—, 

ox ay oz 



(5) 



(G) 



L = xp y - yp x . (7) 

We define also 

P x = Px + By/2, P y =p v -Bx/2. (8) 
Then the operators P x , Py and L each commute with H and p z . But they don't commute with each other: 

[P X ,P V ] = ihB, [L, P x ] = ihPy, [L, P y ] = -ihP x (9) 
We may choose the wave-function tp(x,y,z) to be an eigenfunction of p z with eigenvalue Hk z : 

i>{x,y 1 z) = e lk ^<j ) {x,y). (10) 



HiP = e lk * z [H' + h 2 k 2 z /{2m)](t), 



(11) 



and from now on we will be concerned with H' and (f>(x,y). 

We may choose 4> to be an eigenfunction of any one of the three operators P x ,P y ,L. We discuss the two 
cases (i) eigenfunction of P x (by rotational invariance we can equally well choose it to be an eigenstates of 
P y ) and (ii) eigenfunction of L. These wave-functions may be written in a more compact form if we define a 
magnetic length / by 



In case (i), we have 



0fe x ,n(a;,y) 



'2tt 



f h/B. 



x ixy /2l' 



€n(y), 



where hk x is the eigenvalue of P x , and £„ satisfies 



2m 



-^ + {(y/i 2 )-kxy 



(12) 



(13) 



(14) 



E' being the eigenvalue of H' . This equation describes a simple harmonic oscillator whose equilibrium point 
is shifted, and the eigenfunctions have the form 



where n is a positive integer, 



Uy) = i- 1/2 un[(y-v)/i], 



V = l 2 k x , E' = —(n+h 
m I 



(15) 



(16) 
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and u n (w) is the real, normalized solution of 

d 2 



dw 2 



w 



u n (w) = (2n + l)u n (w). 



(17) 



Next we take the case (ii), where we use energy eigenfunctions which are also eigenfunctions of angular 
momentum L with eigenvalue HM . We denote these normalized wave-functions by 



(18) 



Eigenfunctions of type (i), in equation (|13|) . are non-normalizable and have a continuous degeneracy, whereas 
those of type (ii) in (|18f) are normalized and have a discrete degeneracy. Nevertheless, it should be possible to 
express each type in terms of the other. To this end, we consider a superposition of the form 

(19) 



(20) 



XnJ = j dk x f(k x )<l)k x ,n{x,y) 

where / is some function to be determined and 

dxdyXn' jXnj =<W J dk x \f(k x )\ 2 . 

In order to make (|19|) an eigenstate of L, we require 

hMxnj = Lxnj 



= -ih(2nl)- 1 / 2 J dk x f(k x ) b|-y^j 



x exp[ik x x - ixy/2l 2 ]u n [(y/l) - lk x }. (21) 

We convert x and x 2 into derivatives of e lk " x and integrate by parts to put these derivatives onto / and u. In 
this way, we find that equation (|21|l is satisfied provided / obeys 



I 2 dkl 



+ l 2 k 2 J = [2(n - M) + 1]/. (22) 

This equation is consistent, and its normalised solution is just 

f{k x ) = l x ' 2 u n ^ M [lk x ]. (23) 
Thus, finally, the required eigenfunction of L is determined to be 

Cn,M(z,y) = l 1/2 / dk x u n -M(lk x )4>k x ,n(x,y). (24) 



Using the ortho- normality of the u n coefficients in (|24|l . we can invert this equation to find the momentum 
eigenfunctions in terms of the angular momentum eigenfunctions: 

n 

(t>k^,n(x,y) = l 1/2 ^ U n -M(lkc)Cn,M(x,y). (25) 

A-/=-oo 

The text-book treatment (see for example 2J) of the diamagnetism of a free electron gas is usually formulated 
in terms of the eigenfunctions of momentum 113fl . The sample is considered to have a finite size with dimensions 
S x , S y in the x- and y-directions, in which case states contribute only for 

\V\ < S v . (26) 

So the number of states with a given energy and value of k z is 

BS x S y _ S x S y . 
2irh ~ 2nl 2 [ ' 
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It ought to be possible to carry out the analysis using the angular momentum eigen- functions i|18fl . In that 
case, it is natural to consider a cylindrical sample, axis along the z axis, with radius R. The asymptotic form 
of the wave function (fTH|) . for large p, is proportional to 

e iM V"~ M exp(-p 2 /4Z 2 ), (28) 

and this has its maximum at 

Pmax = v/2(2n - M) I. (29) 

Thus we expect states to contribute for 

2n-M<|j. (30) 



Then, for given n, the range of M is 



2n - |^ < M < n, (31) 



so that the number of possible values of M is 

^ - n. (32) 

The relevant values of n are of order (mkTl 2 / H 2 ) which is small compared to (R 2 /l 2 ) for typical values of B, 
R and T (temperature). Thus l|3*2l gives a factor 

— (33) 

which is proportional to the area just as in (|27|) . 



III. WEAKLY VARYING MAGNETIC FIELD 



We now allow the magnetic field to vary slowly with z, but we restrict ourselves to motion near the z-axis. 
Along the axis, we use an expansion 

B z (z) = B [l + hiz/a) 2 + b 2 (z/a) 4 + ...] (34) 

where a is a characteristic length. 

A realistic case would be a pair of two similar current loops, each perpendicular to the common z-axis and 
of radius a, separated by a distance 2c?. This system, known as the Helmholtz coil, is shown in Fig. 1. 




Figure 1: Schematic illustration of a Helmholtz coil. 



Near the origin the field has an expansion of the above form, where Bq is determined by the current carried 
by the loops and (in terms of r = (d/a) 2 ) 

3(4r-l) 15(1 - 12r + 8r 2 ) 

bl = WTW' 2 = W+W ' (35) 
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Off the axis, the magnetic field must have some radial component. Up to second order (that is the b\ term), 
we take the vector potential to be 

A = is [l + &i(z 2 -p 2 /4)/a 2 ](-y,M). (36) 

This form follows by rotational symmetry (if the field is generated by currents in circles about the axis), 
together with the gauge choice V • A= and Maxwell's equation V 2 A = 0. 

As the leading approximation, we keep just the b\z 2 term in (|34[1 (we shall see that the p 2 term is effectively 
smaller). We seek an energy eigenfunction of the Hamiltonian 0, with B now z-dependent, of the form 

W{z)*(x,y). (37) 



Then the equation for W is (we now define I to be y/h/Bo) 

d 2 W ,„ Mz 2 , TT 2mAE„ T 
-^ + i2n + 1) ^ W =— W > (38) 

where the second term on the left comes from E' in (|14H when B is z-dependent. This is just the equation for 
a simple harmonic oscillator and, for k z = 0, the energy associated with the z motion is 

ft2 i 

AE = -y/bi(2n+l)(N+~). (39) 

mal I 

We can now treat the z 4 term in 134|l and the p 2 term in i|3t)|) as perturbations. These terms are of the same 
order. This is because the transverse size of the wave-function is of order I whereas the longitudinal size is of 
order \/Ta. Note that, as an example, for a magnetic flux density of 1 tesla, I ~ 2.5 x 10~ 8 m. 

Treating the z A term by perturbation theory, we find an additional energy 

aE ' = 4SM2^T)< 1 + 2w + 2Jv2 >' < 40 » 

Note that this is independent of B and usually very small if a is macroscopic. However, for large values of N, 
the contributions from the 62 term and higher terms in l|34|) may be significant. 

We also treat by perturbation theory the p 2 term in 1)36(1. To this end we require, inserting the p 2 term in 
(|36[) into the Hamiltonian (JSJ), 

AE " = - , o m \„ I dxdy^ nM (2Lp 2 + hp'/l 2 )^. (41) 



\Qma 2 l 2 j 

These expectation values can be worked out using: 



As a result, we have 



dxdy^p 2 ^ = (2n - M + (42) 
dxdy^*p A ip = [{2n -M+ I) 2 + 2n(n - M) + 2n-M + l]l 4 . (43) 



h 2 b 

AE" = - \ [6n(n + 1) - M(M + 1) + 2(1 - nM)]. (44) 



This is independent of B and, in general, is much smaller than l|40|) for large values of N. 



IV. NON-PERTURBATIVE APPROACH TO A SLOWLY VARYING FIELD 

We have seen that when the quantum number associated with the z-motion is large, N >> n, the dependence 
of B z upon p may be neglected in first approximation. In this case, the vector potential may be written as: 

A=^B z (z)(-y,x,Q). (45) 
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This leads to the Schrodinger equation: 



d 2 ip I dtp ( M 2 eB z M e 2 B 2 



dp 2 p dp 



he 



Ah 2 c 2 1 



d 2 ^ 



2m 



dz 2 



h 2 



(46) 



where we have taken the energy eigenfunction to be also an eigenfunction of the angular momentum L z with 
eigenvalue HM. Let us consider a solution of the form: 



ip{(t>,P,z) = F{<l>,P,z)f(z) 
and assume, since the field is slowly varying, that: 



dF 




dF 


dz 


« 


dp 



Then, the Schrodinger equation l|46l) takes the form 
~d 2 F 1 dF 



1 
F 



dp 2 p dp 



n 2 



eB z M e 2 B 2 2 



he 



Ah 2 c 2 



Id 2 / 2m 



fdz 2 



h 2 



E 



The right hand side of (|49|l depends on z only, so that we have: 

Itz 2 - + w Ef - -w E {z)f 



and 



d 2 F 1 dF 
dp 2 p dp 



M 2 eB z M e 2 B 2 p 2 
1 f 1 £i_ 

n2 



he 



Ah 2 c 2 



F 



2m 



E'(z)F. 



But we know how to solve (|5ip. in which z is just a parameter. Comparing with we see that: 



E'{z) 



eh 



:)B,(z) 



Thus, H5Uf) may be written as: 



d 2 f 2m 
dz 1 + Ji 2 ~ 



E- 



eh(2n + 1) 
2mc 



B z (z) 



f = 0, 



(47) 



(48) 



(49) 



(50) 



(51) 



(52) 



(53) 



which is just an ordinary differential equation. 

We can now check the assumption l]48fl. The solutions of l|51() have the form l|18|). namely: 



F(<j),p,z) = ( n ,M( 



eB z 
2hc 



P) 



Hence, the condition l|48() requires that: 



p dB z 



« 1. 



Although this cannot be satisfied for all p, we note that the solution Ij54(l contains the factor: 

r eB * 2i 



so that typical values of p are of order ^JhcjeB z . Thus, we expect l)55|l to be valid for: 

1 dB z 



e B° z 



3 / 2 



« 1 



(54) 



(55) 



(56) 



(57) 



This can be satisfied provided we take B z to be varying slowly enough. For instance, if B z (z) has the form 
(|3"P with \z\ of order a, then the left hand side of ijSTJl is of order l/a, which is very small for typical values of 
Bo and a. 
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V. QUANTIZED ENERGY LEVELS IN A HELMHOLTZ COIL 

A simple example of a magnetic mirror is provided by the pair of current loops shown in Fig. 1. Then, as 
can be seen from equations (|34|) and l|35fl for 2d > a, the magnetic field can increase enough in the region near 
the current loops, so that a charged particle may eventually be reflected out of this region towards the center 
of the coil. The classical restoring force which may confine the particle along the z-axis is provided by the 
p-component of the magnetic field: 



B(*) 



1 dB z 



2 9 ^ 



B 7 i. 



(58) 



which can be derived from (|45|l using B = V x A. 

In this section, we will discuss the quantized energy levels of a charged particle in a Helmholtz coil. To this 
end, we will analyse the basic equation l|53|l . where 



B z (z) = ]-B Q {d 2 + a 2 ) 3/2 \- — ^ ^777 + 77 777 77777 

2 v ' {[(z + d) 2 + a 2 } 3 ' 2 [{z - d) 2 + a 2 } 3 ' 2 

The general behavior of this field, for 2d > a, is shown in Fig. 2. 



(59) 




Figure 2: Pattern of the field B z (z) in a Helmholtz coil. Circles indicate the classical turning points. 

Since (|53|) was derived under the assumption that N >> 1, our analysis will be more appropriate in the 
semiclassical regime. This may be studied conveniently in the WKB approximation, which is particularly useful 
since we are dealing with a slowly varying potential: 

V(z) = ^(2n+l)B z (z). (60) 

One of the most interesting results arising from the WKB method is a semiclassical estimate for the quantized 
energy levels in a potential. Matching the WKB wave function at each of the classical turning points, which 
are determined by the relation E — V(z), leads to the Bohr-Sommerfeld quantization condition: 

1 



y/2m[E - V{z)]dz = (N+ -)irh , (61) 
z + 2 

where the classical turning points for bound motion, ±Z+, are situated inside the well as shown in Fig. 2. (We 
neglect, in first approximation, the very small probability of tunneling through the potential barrier.) 

In general, equation (|61|l is rather complicated and can only be solved numerically. However, when d/a is 
somewhat larger than 1/2, it may be solved in closed form, since in this case z/a will be effectively of order 
1/2 or smaller. Then, we can expand the potential V(z) [see (|59|) and (|6L)fl ] up to terms which are of quartic 
order in z/a, as shown in equations (|34[1 and (|35|) . In this approximation, the quantization condition (|61|l may 
be written in the form: 



_z+ V a a V (2^ + l)ei>o 2 



where the dimensionless parameter £ is defined by: 



£ = (2„ + 2 Tw £ - < 63 > 
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We can now determine explicity the positions of the classical turning points, which are given by: 

-6i± v/6?+4b 2 (£-l)i 1/2 



Z± = a 



2b 2 



(64) 



where ±Z_ are the turning points for unbound motion, which are situated outside the well as shown in Fig. 2. 
Then, the integral appearing in (|62jl can be evaluated in closed form in terms of the complete elliptic integrals 
of the first (F) and second (E) kind [3j, with the result: 



3a 2 



1 ch l 

(2n+l)e_B + 2 



(65) 



This is a transcendental equation which determines implicitly the energy in terms of the quantum number N. 
The energies of the bound states have an upper bound E max which corresponds to the maximum value of the 
potential. At this energy, Z + = Z_, so that using l|63[l and (|64|l we obtain: 



E n =-Z O 1 

Zmc 



1)1- 



Abo 



(66) 



When Z + = Z-, the equation l|65|l simplifies considerably and fixes the maximum value of the quantum 
number N , which is given by the relation: 



N, 



1 



y /2(2n+l)b 3 1 /2 a 
Sir \b 2 \ I 



max H~ 2 



(67) 



where I — \J ch/eBo is the magnetic length. Hence, N max is in general a very large number, which is in 
acordance with our previous assumption. 

A more explicit relation giving E as a function of N can be obtained by solving numerically the equation (|65|) . 
As an example, the numerical solution for n — and d/a — 3/5 is shown in Fig. 3. Here, the numerical values 
of £ max = 1.028 and N max — 1.12 x 10 6 are in good agreement with the corresponding results obtained from 
the closed form expressions (|5T)ll and l(57)l. One can see that the numerical solution can be fitted reasonably 
well, for large quantum numbers, by a simple phenomenological form like N 5 / & . 




Figure 3: Numerical solution (crossed) of Eq. I|65|l for the energy spectrum and plot (solid) of the form iV 5 / 6 for large 
TV. 

Finally, we note from and (J35J that for d/a = 1/2, bi = and b 2 = -144/125. Then, the system of the 
two current loops provides a practically uniform field in the central region of the Helmholtz coil. In this case, 
H66J) reduces to the well known expression for the Landau levels which occur in a constant field. Therefore, 
one may regard the quantized energies described by (|65(l and 166(1 , as an extension of the Landau levels to the 
case of a slowly varying magnetic field. 



VI. CONCLUSION 



In this paper, we have discussed several aspects concerning the quantum behavior of a charged particle in 
a static magnetic field. We have treated the issue of the relation between the choice of gauge and the choice 
of the diagonal operators which commute with the Hamiltonian of the system. We have also developed some 
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approaches which may be useful for physical applications in a slowly varying magnetic field. These methods 
have been applied to study the quantized energy levels of a charged particle in a weakly varying field. Such 
a field may be present, for example, in a magnetic mirror like the Helmholtz coil. We have shown that this 
energy spectrum represents an interesting extension of the well known Landau levels which occur in an uniform 
magnetic field. 

This work was supported in part by US DOE Grant number DE-FG 02-91ER40685, by CAPES, CNPq and 
FAPESP, Brazil. 
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